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Abstract. We propose a quasiparticle approach allowing to compute the equation of state of a generic 
gauge theory with gauge group SU(A^c) and quarks in an arbitrary representation. Our formalism relies 
on the thermal quasiparticle masses (quarks and gluons) computed from perturbative techniques, in which 
the standard two-loop running coupling constant is used. Our model is minimal in the sense that we do 
not allow any extra ansatz concerning the temperature-dependence of the running coupling. We first show 
that it is able to reproduce the most recent equations of state computed on the lattice for temperatures 
typically higher than 2 Tc- Well above Tc indeed, an ideal gas framework with thermal masses is indeed 
expected to be relevant. Then we study the accuracy of various inequivalent large- A^c limits concerning 
the description of the QCD results, as well as the equivalence between the QCDas limit and the J\f — 1 
SUSY Yang-Mills theory 
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1 Introduction 

More than two decades after pioneering works [T| , the phe- 
nomenology related to the deconfined phase of QCD, i.e. 
the celebrated quark-gluon plasma (QGP), is still a mat- 
ter of intense investigations, both experimentally and the- 
oretically. On the experimental side, the QCD matter is 
or will be studied in heavy-ion collisions at RHIC, SPS, 
FAIR, and the LHC. We refer the interested reader to the 
reviews [2[ for a survey of the main experimental issues in 
QGP physics. 

On the theoretical side, the study of gauge theories 
at finite temperatures, among which QCD plays a central 
role, deserves much interest since it is both a challeng- 
ing problem in itself and possibly applicable to current 
experimental results. In the present study, we will focus 
on the two following approaches: lattice QCD and quasi- 
particle models. The lattice formulation of QCD allows to 
extract nonperturbative informations from QCD through 
numerical computations, and lattice data are often used 
to fit phenomenological model's parameters when no ex- 
perimental result is available. An observable that can be 
accurately computed on the lattice and by resorting to 
effective approaches is the QCD equation of state (EoS), 
i.e. the evolution of various thermodynamical quantities 
versus the temperature. Although not directly measurable 
from experiments, the EoS is of crucial importance to un- 
derstand the phase diagram of QCD. One can for example 



read the reviews [3] for complete discussions of the vari- 
ous EoS computed so far: with or without quark flavors, 
at zero chemical potential or not, etc. More results can be 
found for example in [31[S]. 

Starting from a more phenomenological point of view, 
quasiparticle models are based on the assumption that the 
QGP can be seen as a gas of deconfined quarks and glu- 
onfQ. Pioneering works [71[51ll] have shown that the con- 
dition for quasiparticle models to reproduce the various 
EoS computed in lattice QCD above the critical temper- 
ature Tc is that the quasiparticles gain thermal masses. 
Those ad hoc masses had to be fitted on the lattice data. 
More recently [TUlfTT] . Hard- Thermal-Loop (HTL) tech- 
niques have been used to compute analytically the behav- 
ior of thermal masses at leading and next-to-leading order 
in the coupling g from QCD. As it will be recalled in Sec.[5J 
such HTL formulas are of interest to build a more firmly- 
based quasiparticle formalism, see e.g. [12] . Our aim is to 
use a formahsm as simple and general as possible so that 
we can compare it to a maximal number of recent lattice 
results. This will be done in Sec. [3] for a pure Yang-Mills 
plasma and in Sec. H] for a general QGP, with light (and 
heavy) quark flavors, and a nonzero chemical potential. 
Moreover, we will not explicitly particularize our formu- 
las to QCD as usually done: Taking into account a SU(A'c) 
gauge group and quarks in an arbitrary representation of 
that group, we will be able to study various large- TVc lim- 
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^ Remark that a quasiparticle description of hypothetical 
quark stars has been proposed historically as an application 



of the quark hypothesis, see [6]. 
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its of our model and to check their accuracy to reproduce 
the QCD results. To our knowledge, that problem has not 
been addressed yet in quasiparticle models and is stud- 
ied in Sec. [SI as well as the case of the — 1 SUSY 
Yang-Mills EoS, that will also be computed in that sec- 
tion. Notice that the question of knowing "how far TVc = 3 
is from oo" is still widely studied nowadays (see e.g. [13'), 
and has gained a renewed interest, in particular thanks to 
AdS/CFT-inspired approaches of QCD; one can quote for 
example the review jl4J . After some comments on heavy 
mesons, we summarize our results in Sec. [HI 



2 Quasiparticle model 

Let us begin by recalling the results that will be used in 
the present study. It is worth saying that those formulas 
are not new in themselves: they have been obtained in 
previous studies by using perturbative methods. The pe- 
culiarity of the present work is rather to apply them with- 
out introducing phenomenological modifications, with one 
set of parameters for all the possible cases considered. 



2.1 Analytic results 

First of all, the two-loop running constant reads [TS] 
87r2 /3i In [2 In (5)] 1 



1 



/3oln(i)l ln(i) 



(1) 



where the temperature, T, has been chosen as the typical 
energy scale of the system, while A (or Aqcd) is the free 
parameter of the model. f3o and f3i are the coefficients of 
the two-loop /3-function, reading, for a theory with gauge 
group SU(A'^c) and nj light quark flavors in the represen- 
tation R of that group. 
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In the above equations, Cu is the value of the quadratic 
Casimir operator in the color representation R, and T^^ = 
C}idimii/{N^ — 1), where dim^ is the dimension of the 
representation R. A general method to compute Cr can 
be found for example in |16) . 

The crudest way of modeling a QGP would be to as- 
sume that it is an ideal gas made of gluons (Bose-Einstein 
distribution) and of Nf quark flavors (Fermi-Dirac distri- 
bution), all assumed to be massless. Therefore, the dis- 
persion relation is e{k) — \k\. In that case, assuming that 
each quark flavor is in the same color representation R and 
has a nonzero chemical potential fi, the pressure would be 
given by the Stefan-Boltzmann-type law 
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T)+Li4(-eT)] dimflA^/}, (3) 



where liis{x) is the poly logarithm function of order s. 
However, if ^ was the QGP pressure p, the ratio p/T^ 
should be constant in disagreement with all the lattice re- 
sults concerning the QGP EoS. Strictly speaking, it has 
been shown that this last ratio is lower than the Stefan- 
Boltzmann case at any T < 5Tc, while reaching p « 
PSB would demand huge and probably phenomenologi- 
cally irrelevant temperatures: around (10^, 10^) x Tc fol- 
lowing |17| . The necessity of going beyond the Stefann- 
Boltzmann limit has appeared soon after the first lattice 
computation of the EoS in pure glue SU(2) theory [T^. In 
the phenomenological studies [71IH] , it has been shown that 
an agreement between the lattice EoS and the one com- 
puted assuming that QCD above Tc consists in an ideal gas 
of quasiparticles (quarks and gluons) can be reached only 
if the quark and gluon masses are temperature-dependent. 
Those thermal masses have to be fitted on the lattice data 
if one wants an agreement in the whole temperature range 
studied by lattice calculations, typically (1.0,4.0) x Tc. 

Instead of introducing ad hoc thermal masses, it is rel- 
evant to recall that perturbative techniques, or HTL re- 
summation methods, allow to compute the thermal quark 
and gluon masses at large temperatures from diagrams of 
self-energy and vacuum polarisation-type, see for exam- 
ple [T?[f201[^[^ . It appears from those studies that the 
thermal gluon and (anti) quark square masses are given by 



V 9^ rTi2 



{Nc + TRnf)2-r' + ti^^TRnf 



rj-iZ 



(4) 



where a nonzero chemical potential fi is associated to the 
light quark flavors. Notice that in the following, Nf will 
denote the total number of quark flavors, while Uf is the 
number of light quarks. The Tr factor is often absent from 
the formulas found in the literature since it is simply equal 
to 1/2 for fermions in the fundamental representation. 

The formulas ^ are valid at the leading order in 
g, thus one has to be aware that our model may only 
be applied at temperatures well above Tc. An estima- 
tion of next-to-leading order corrections using HTL tech- 
niques has been made in 22 : The corrective terms 6nri = 



'Ncg'^Tmg/V2Tr and Or 



-Cr g^TiTLg/ \/27r have been 



found. For the sake of simplicity, we nevertheless choose 
to keep the formulas Q in our formalism, as done in 
most of the quasiparticle studies. A possible justification 
is that such corrections mainly weakens the increase of the 
leading-order thermal masses with T. Provided that the 
model is not applied within a too large temperature range, 
as we will do, that feature can simply be mimicked in first 
approximation by a larger value of A, our free parameter. 



2.2 Thermodynamical properties 

The grand-canonical partition function Z of an ideal quan- 
tum gas containing quarks and gluons whose thermal masses 
are given by (|1]) can be obtained using standard statistical 
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mechanics. The pressures associated to each species can be 
computed thanks to the standard relation p — TlnZ/V, 
V being the volume of the system. In units where h — c — 
ks — 1, the gluon pressure reads (see e.g. [13]) 



(5) 

while the contribution to the pressure of one quark flavor 
with chemical potential /i is given by 



T 

Pq{mq) = 2dimfl — 



dk fc2 In 1 



(6) 

the only change occuring at the level of the chemical po- 
tential. The degeneracy factor in the above expression in- 
cludes the spin and color degrees of freedom. We remark 
that a possibly nonzero bare quark mass, rfig, has been 
quadratically added to the thermal quark mass - thermal 
effects come indeed as a modification of the bare quark 
propagator. The corresponding antiquark pressure is given 

by 



Pq{mq) = 2 dimn / dkk'^\n\l + e 



(7) 

In general, our results will be normalized to the Stefan- 
Boltzmann pressure ^ so that they remain finite at large 
Nc- Notice that psB given by Q is nothing but the sum 
Pg + Nf(pq + pq) in which rUg = iriq = fhq ~ 0. 

In the following we will mostly focus on the trace 
anomaly, which is defined as A — e ~ 3p and is a par- 
ticularly relevant observable since it is a direct measure of 
the nonideal character of the deconfined medium. It can 
be computed from the pressure by using the thermody- 
namical relation A = T^drip/T^), or equivalently 



PSB 



\PSB J 



(8) 



It has been previously shown that such a quasiparticle 
picture is able to reproduce some lattice data from Tc to 
higher temperatures, provided that the T-dependence of 

is modified, typically by replacing ln(r/A) in g^ by 
some function ln(/i(T/yl)), with h{x) oc x at large x, see 
for example [9l[12]. In the present work, we will adopt a 
different point of view. We think indeed that it is not 
worth modifying the above perturbative results in order to 
obtain a model in agreement with lattice QCD in the range 
(1.0, 2.0) X Tc. The reason is that, at those temperatures, 
the QGP is far from being an ideal gas, as shown by the 
maximal value of A/T'^, systematically reached around 
1.2 Tc in lattice computations. This means that, in that 
range, not only the behavior of should be modified, but 
mostly the framework itself, which needs to go beyond the 
ideal-gas assumption. That is why we will apply the above 
formalism as it is, but where it is the most likely to be 
trusted, that is at temperatures well above Tc, typically 
above 2Tc. The arbitrary is then considerably reduced. 



3 The SU(A^c) gluon plasma 

The simplest case to which the present approach can be 
applied is the gluon plasma, whose EoS has been first com- 
puted in SU(2) il8j and then in SU(3) lattice QCD in g4], 
while more accurate results have recently been obtained 
in T using gauge groups ranging from SU(3) to SU(8). 

In the pure glue case, Nf = rif = and p — Pg. Then, 
it can be observed from Eq. ^ that the thermal gluon 
mass is independent of Nc since nig oc Ncg"^. This last 
term is constant with respect to Nc since g^ — X/Nc from 
(HI; A is actually the 't Hooft coupling. Consequently, pg 
is schematically equal to {Nc — 1) f{T), leading to ratios 
p/psB and A/psB which are independent of A^c- That fea- 
ture of the gluon plasma has indeed been observed in [1] , 
where the EoS for the different gauge groups considered 
are found to be compatible up to the error bars (see also 
Fig. [T|). Notice that the above discussion is based on the 
leading-order formulation oimg. In principle, higher-order 
terms might introduce other dependences in Nc] typically 
one could expect nig to be constant up to terms in 1 /Nc as 
it is generally the case with gluonic observables. Such sub- 
leading corrections cannot be studied within our frame- 
work, where we focus on the dominant behavior of the 
EoS. We thus stress that, in the rest of this paper, the 
large- A^c behavior of the thermodynamical quantities is 
meant to be that at dominant order in Nc, i.e. our con- 
clusions are valid up to subleading corrections in A'^c- 

The general framework presented in the previous sec- 
tion, when particularized to the present case, describes an 
ideal gluon gas, and thus may be applied at temperatures 
where the interactions between gluons are not dominant 
with respect to their kinetic energy. Looking at the lattice 
results of [31 in Fig. [U it appears that the gluon plasma is 
maximally far from an ideal gas around 1.2 Tc. Then A/T'^ 
quickly decreases, but our approach should clearly not be 
applied at temperatures lower than typically 2 Tc, where 
A/T^ becomes well smaller than its maximal value. The 
coherence of that discussion would require the interaction 
potential between two gluons above 2 Tc to be signifi- 
cantly smaller than around Tc. Informations about that 
potential can be deduced from lattice results concerning 
the free energy between static quarks: The free energy be- 
tween two gluons should indeed be the same as the one 
between a quark-antiquark pair up to an overall color fac- 
tor in virtue of the Casimir scaling [25] . The gluons being 
light particles, they should be especially sensitive to the 
large-distance part of the interaction potential, and it ap- 
pears from lattice calculations that the asymptotic value 
of the free energy between a quark and an antiquark at 
2.0 Tc is less than two times the one at 1.2 Tc [IHj- So it 
seems reasonable to trust our minimal quasiparticle model 
mostly well above Tc, where the color interactions are sig- 
nificantly weakened. 

The only parameter having to be fitted is the energy 
scale A, present in g^. That scale is usually interpreted as a 
typical energy scale for the confinement. Since the temper- 
ature T is used as an average energy for the system under 
study, it can be guessed that A/Tc — 0{1). Indeed, the 
deconfinement can intuitively be expected to occur when 
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Fig. 1. (Color online) Pure glue trace anomaly, normalized to 
the Stefan-Boltzmann pressure ^ and computed by using ((SJ 
with p = Pg given by ^ (solid line). The result, independent 
of Nc and computed with A/Tc — 0.88, is compared to the 
corresponding lattice data taken from ^ for the gauge groups 
SU(3) to SU(8) (full symbols), and from [18] for SU(2) (empty 
symbols). 



the average energy of the thermal fluctuations is able to 
"compete" with that of the confinement, thus when T be- 
comes similar to A. Here A is fitted to the SU(3) trace 
anomaly computed in [4] in the range (2.0,3.5) x Tc, that 
is the case for which a maximal number of data is avail- 
able. In accordance with our guess, we obtain A = 0.88 Tc, 
a value that leads to an excellent agreement with the lat- 
tice data as shown in Fig. [TJ We have plotted our model's 
results up to 1.5 Tc also, and it appears that the model 
begins to miss the data at such low temperatures, as ex- 
pected. Notice that, if the lattice data are compatible with 
our model for all Nc at high temperatures, the SU(2) lat- 
tice data are clearly distinct from the other values of Nc 
in the range (1.0,2.0) x Tc- Two interpretations of that 
fact can be done, not excluding each other. First, those 
data, coming from |18| . are more than twenty years old. 
The lattice used was clearly smaller than the most recent 
calculations [4] , and it is possible that what we observe is 
actually a lattice artefact. Second, the SU(2) data begin 
to differ significantly when we have argued that neglect- 
ing the interactions between gluons becomes irrelevant. 
Those interactions should induce corrections, probably in 
inverse powers of l/iV^ as it is usually the case for gluonic 
observables, that can be significant when Nc — 2. Only a 
new lattice calculation of the SU(2) EoS could clarify that 
issue. 

An extension of the present framework to temperatures 
near Tc would demand the inclusion of gluon-gluon inter- 
actions within a non-ideal Bose gas framework. A first step 
in that direction can be found for example in p71f25] . How- 
ever, such a question will not be addressed in the present 
work, which is an attempt to find a model describing most 
of the currently known lattice data with as few parameters 



as possible, namely A, Tc and the quark masses. Concern- 
ing the interactions nevertheless, it is worth mentioning 
that two gluons can be found in the following color chan- 
nels: 



Nc-1 



□ □ 




where all the Young diagrams of the first (second) line 
have Nc — 1 {Nc — 2) rows. Notice that the last Young dia- 
gram is present only if A^c > 3, and that only the first three 
diagrams exist at A'c = 2. Any corresponding potential of 
one-gluon-exchange form, typically a Yukawa potential, 
would be proportional to the color factor {Cr^ —2Nc)g'^/2, 
Rg being the color representation of the gluon pair. More- 
over, it should involve the Debye mass, linked to the ther- 
mal gluon mass. Recalling that = X/Nc, the color fac- 
tors corresponding to the Young diagrams ^ read (di- 
vided by A) -1, -1/2, 1/Nc, 0, -1/2, 0, and -1/Nc re- 
spectively. Since the thermal gluon mass is constant with 
Nc at leading order, this means that the dominant gluon- 
gluon interactions are independent of A'c- So, from a quasi- 
particle picture, one expects the pure glue EoS to be con- 
stant with Nc in the whole range T >Tc at the dominant 
order, as observed first in [3] and confirmed in |29| . 

We mention for completeness that the thermodynami- 
cal properties of pure gauge QCD has been recently com- 
puted in HTL perturbation theory at three loops [30] ; the 
results compare very favorably with lattice QCD down to 
(2, 3) X Tc and the independence of the results with respect 
to Nc is readily checked. Also, nonperturbative effects like 
a contribution coming from the gluon condensate could 
be taken into account, as suggested for example in [3n[5^ 



4 Adding quark flavors 

4.1 Vanishing chemical potential 

So far we have shown that our simple quasiparticle ap- 
proach is able to reproduce the lattice EoS of a gluon 
plasma well above Tc for any gauge group provided that 
the parameter A is fitted to the data. The next step is 
now to add the quark flavors to our formalism and see 
whether the corresponding lattice data can be reproduced 
or not in the same temperature range. We choose to keep 
the value of A previously fitted to the pure glue case so 
that the arbitrary of the model is reduced, while the quark 
masses are taken from the PDG [M] . The only needed ex- 
tra parameter is the QCD critical temperature. Lattice 
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simulations generally find to be located in the interval 
(150, 200) MeV when light quark flavors are present [35l 
[5M57] . Here we take a quite typical value of 180 MeV, but 
our results would not be significantly affected using any 
Tc located in the above interval. The parameters that will 
be used in the rest of this paper are 

yl = 0.88Tc, 180 MeV, 

TO„ = mrf = OMeV, = 150 MeV, 
fhc = 1270 MeV, fhb = 4200 MeV, (10) 

where the u and d quark masses, given by 2.5 and 5 MeV 
respectively in the PDG, have been taken equal to zero. 
In the lattice calculations [51 1551155] , to which we will com- 
pare our results, the strange quark mass has nearly its 
physical value, rhg ~ fhs (the hatted symbols denote the 
lattice masses), while the u and d masses are such that 
rhu = rhd = mg/lO, that is a larger value than the PDG 
ones. Moreover, in [51l40j . the heavy quark masses are such 
that rhc/rhs = 10 and rhb/rhs — 38, that is a charm mass 
compatible with the experimental one, while the lattice 
bottom mass is quite heavy, around 5700 MeV. We have 
checked that the curves we obtain are nearly indistinguish- 
able if the u and d quark masses are allowed to vary from 
to 20 MeV, i.e. values close to the experimental and 
lattice ones respectively. 

Most of the current effort in lattice QCD has logically 
been devoted to the study of a SU(3) QGP with Nj — 
Uf = 2-\-l light quark flavors and zero chemical potential, 
whose pressure reads in our framework 

Plq =Pg + 2pg(0) +2pq{0) + Pq{fhs) +Pq{fhs), (11) 

with Nc — 3 and Nf — Uf = 2 + 1. The corresponding 
trace anomaly, computed thanks to Eq. ([S]) in which piq is 
used, compares very well, when normalized to psB, with 
the available lattice data from the most recent computa- 
tions, involving lattices with a larger number of tempo- 
ral slices Nt = 6 and 8 (see Fig. ^ . Notice the following 
numerical coincidence: When they are expressed in units 
of Tc and normalized to the Stefan-Boltzmann pressure, 
those data are very close to the ones of the pure glue case. 

Because of the statistical suppression in e"™/-^ for a 
given particle species [see Eq. (jG])], the dominant contri- 
bution to the EoS is expected to come from the lightest 
quarks. However, at high enough temperatures, the ra- 
tio m/T can become small even with heavy quarks, al- 
though the decrease is slower than for light quarks. That 
is why the contribution of deconfined heavy quarks should 
be taken into account in more realistic approaches of the 
quark-gluon plasma at high temperatures. The charm and 
bottom quark contribution to the EoS has, to our knowl- 
edge, only been investigated on the lattice in [51HD]. These 
recent results can be compared to our model's prediction, 
reading 

Phq = Plq +Pq{mc) + Pq{mc) 

for iV/ = 2 + 1 + 1, 
= Plq +Pq{rnc) +Pq{mc) +Pq{rhb) +Pqifhb) 

for TV/ = 2 + 1 + 1 + 1, (12) 
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Fig. 2. (Color online) Trace anomaly of a SU(3) quark-gluon 
plasma with 2+1 light quark flavors at zero chemical po- 
tential, normalized to the Stefan-Boltzmann pressure Q and 
computed by using (O with p — piq given by (solid line). 
Parameters pO[) are used. The result is compared to the cor- 
responding lattice data taken from [38] for At = 4, 6 (empty 
symbols), and from [39] for At = 6, 8 (full symbols). The data 
of both works are indistinguishable for At = 6. Pure gauge 
results for the group SU(3), taken from [5, are plotted for 
comparison (crosses). 



with Nc — 3 and Uf = 2 + 1 in both cases. The maximal 
considered temperature is 5 Tc « 900 MeV: The heavy 
quark loops can be neglected at such a typical energy, 
justifying our choice for nj. As shown in Fig. [3] no lat- 
tice data are currently known above 2 Tc but our results 
seem to extrapolate coherently the available ones. Remark 
that the heavy quark contribution significantly modifies 
the asymptotic behavior of the results obtained with only 
light quarks, as illustrated in Fig. [3] Notice that using 
the lattice b mass of 5700 MeV instead of 4200 MeV only 
slightly modifies the curve obtained at A^t = 2 + 1 + 1 + 1. 
We have not plotted it in Fig. [3] for the sake of clarity. 

4.2 Nonzero chemical potential 

In view of an eventual comparison of quasiparticle mod- 
els to experimental results, the use of a nonzero chemical 
potential is particularly relevant. The most recent lattice 
calculations deal with a nonzero constant ratio fi/T and 
2+1 light quark fiavors .5..41,.42 . Since the error bars 
concerning the change in trace anomaly due to a nonzero 
chemical potential are quite large in these works, we prefer 
to focus on an other observable for which the uncertainty 
in smaller: the change in pressure, defined in our model as 

Notice that, when /i ^ 0, the Stefan-Boltzmann pres- 
sure ([3]) also changes by an amount 

^PSB = Psslf, ~ P5s|^=o • (14) 
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Fig. 3. Trace anomaly of a SU(3) quark-gluon plasma with 
2 + 1 light quark flavors plus charm (solid line) and bottom 
quarks (dashed line) at zero chemical potential, normalized to 
the Stefan-Boltzmann pressure ^ and computed by using (O 
with p = phq given by H12|) . The trace anomaly with only 2 + 
1 light quark flavors is plotted for comparison (dotted line). 
Parameters (|10p are used. The results are compared to the 
corresponding lattice data taken from [5] (symbols). 
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Fig. 4. Change in pressure of a quark-gluon plasma with 2 + 1 
light quark flavors due to a nonzero chemical potential, nor- 
malized to the change in Stefan-Boltzmann pressure (|14|l and 
computed by using ()13|) with p = piq given by (|ll|l (solid lines). 
Parameters (|10p are used. The results are compared to the cor- 
responding lattice data taken from 5 (symbols). Notice that 
the curves obtained at different fi/T are indistinguishable. 



In a previous lattice study [13] , where the EoS is stud- 
ied for 2 + 1 light quark flavors {Nt = 4) at nonzero con- 
stant between 100 and 530 MeV. Among other results, 
the observation is made that the ratio Ap/Apss is con- 
stant, up to the error bars, for any chemical potential in 
the studied range. Following our model, those ratios differ 
by 3% at 2 Tc provided that /i is taken between 100 and 
530 MeV; that difference moreover quickly tends to zero 
as T increases. This confirms that, as first pointed out 
in [43], the dependence in ^ of Ap and ApsB factorizes 
in the same way, thus independently of the quark masses. 
As for the case where heavy quarks are added, only a few 
data are currently known beyond 2 Tc- However, a glance 
at Fig. m shows that the ratios Ap/ApsB, obtained from 
the lattice results of [5], are indeed independent of n/T 
(at Nt fixed), in agreement with our results. Furthermore, 
the present approach seems to agree with the asymptotic 
behavior of Ap/ApsB suggested by lattice QCD. We re- 
mark that the computed curves are logically incompatible 
with the lattice data near Tc, that is where we think that 
our approach can certainly not be trusted anymore. 

5 Large Nc limits 
5.1 The different cases 

Despite its simplicity, the present model is defined for an 
arbitrary gauge group SU(A'^c) and quark representation 
R. Consequently, it allows to check the accuracy of var- 
ious large- A^c hmits in reproducing the EoS in the QCD 
case, that will denote hereafter an SU(3) theory with 2 + 1 
light quarks in the fundamental representation. The only 



criterion to formulate a laige-Nc generalization of QCD 
is obviously that it should reduce to usual QCD when 
Nc = 3; then results can be extrapolated to an infinite 
number of colors and the "closeness" of that limit with 
respect to Nc = 3 can be studied. We mainly find three 
different cases of interest: 

— 't Hooft limit: In that limit, which has been first pro- 
posed in the literature [44 , 45 , the number of colors can 
become arbitrarily large while the 't Hooft coupling 
A = g'^Nc remains constant for any Nc- The quarks 
are still in the fundamental representation, with a fi- 
nite value for A^^. Thus, in that limit, the gluon con- 
tribution to the EoS scales as A'^ while the quark one 
scales as Nc and can be neglected at the dominant 
order. 

— Veneziano limit: That case is similar to the 't Hooft 
limit, excepted that the number of quark flavors is al- 
lowed to become infinite, i.e. Nf = 0{Nc) [H]. This 
implies that the internal quark loops in Feynman di- 
agrams are no longer suppressed at large Nc. In our 
case, the consequence of that limit is that the quark 
contribution to the EoS scales as N"^ instead of Nc. 
We will take explicitly Nf — rif — Nc, which leads to 
three light quark flavors as required for Nc = 3. 

— QCDas hmit: When Nc — 3, the fundamental and con- 
jugate representations are isomorphic. One can then 
in principle define a quark to be in the representa- 
tion i? = 2 = ; this is just a redefinition of the 
terms "color" and "anticolor". At large Nc however, 
the fundamental and two-index antisymmetric repre- 
sentations are no longer equivalent. We call QCDas the 
limit in which Nf remains finite, but with quarks in 
the A2-representation. In large-A'c QCDas the quark 



Fabien Buisseret, Gwendolyn Lacroix: A minimal quasiparticle approach for the QGP and its large- A'^c limits 



7 



loops are not suppressed. Moreover, this case has been 
shown to be equivalent to a large- iVc gauge theory with 
quarks in the adjoint representation |47| . This leads 
in particular to an appealing duality between QCDas 
with iV/ = n/ = 2 and TV = 1 SUSY Yang- Mills; we 
will come back to that point later. 

We mention for completeness the Corrigan-Ramond 
limit [3H]j which has been the first attempt to include 
quarks in higher representations. In that limit, Nf re- 
mains finite and the quark flavors are separated into two 
color representations, either the fundamental one or the 
{Nc — 2)-index antisymmetric one, which are both equal 
at Nc = 3. This allows in particular to represent baryons 
as three-quark state for any Nc, but this limit will not be 
considered in the present paper, where we have chosen to 
keep all the quarks in the same color representation. 



5.2 EoS at large Nc 

We are now in position to study the behavior of the EoS in 
the different large- A^c hmits that we have proposed. The 
results obtained for the trace anomaly of a QGP with only 
light quark flavors are given in Fig. [51 Several observations 
can be made. First of all, the 't Hooft and Veneziano lim- 
its are both excellent approximations of the QCD case, 
and reproduce the lattice data accurately as well. Then, 
the QCDas hmit with Nf — nf — 2 + 1 fails to reproduce 
QCD, and has even a very different qualitative behavior. 
The reason can be understood by studying the general 
structure of the underlying theory, as done in [IS] . Follow- 
ing the conclusions of this last work, based on the study of 
the two-loop /3-function, the gauge theory underlying the 
QCDas hmit is conformal for || < n/ < -^-j and looses 



its asymptotic freedom beyond that value. But, already for 
n/ > ^ = 2.125, it exhibits a Banks-Zaks infrared fixed 
point. Such a feature, not shared by QCD, is thus present 
when Nf = n/ = 2 -|- 1; we think that it explains why 
the physical behavior of the normalized trace anomaly is 
so different in that case than for QCD. If Nf = n/ = 2 
however, no fixed point is present, and it can be observed 
in Fig. [S] that, although being in poor quantitative agree- 
ment, the qualitative behavior of the trace anomaly in 
that limit is at least in agreement with the QCD case. It 
is worth saying also that a Banks-Zaks infrared fixed point 
appears in the Veneziano limit if nf /Nc > 34/13 [IP]; this 
is not the case here since we have taken Nf — nf — Nc- 

It is also of interest to compute the change in pressure 
due to a nonzero chemical potential in the different limits 
of interest: Results are plotted in Fig. [HI As expected from 
the previous discussion, the QCDas hmits are inaccurate 
with respect to the QCD case. However, the differences 
between the 't Hooft and Veneziano limits are more vis- 
ible than for the trace anomaly; the hierarchy of the re- 
sults are coherent with the relative contributions from the 
quark loops in those limits. Consequently, the study of the 
EoS at nonzero chemical potential might be relevant also 
in checking the relative accuracy of different inequivalent 
large- A'c limits of QCD. 



- 't Ilooft 
Veneziano 

- QCDjj,,A?-=n^=2+l 



\ SUSY Yang-Mills 
p4 6, 8 




N =3,^=77 =2+1 



Fig. 5. Trace anomaly of a quark-gluon plasma with light 
quark flavors and zero chemical potential, normalized to the 
Stefan-Boltzmann pressure ([S]), in the different large- Ac limits 
defined in Sec. 15.11 (lines). Parameters (|10|) are used. Lattice 
data taken from |39 , corresponding to a SU(3) quark-gluon 
plasma with 2-1-1 light quark flavors, have been added for 
comparison (circles). 
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W -3,W^-n^-2+l 
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Fig. 6. Change in pressure of a quark-gluon plasma with 2 + 1 
light quark flavors due to a nonzero chemical potential, nor- 
malized to the change in Stefan-Boltzmann pressure (|14p . in 
the different large- Ac limits deflned in Sec. [S](lines). Parame- 
ters (|10|) are used. Lattice data taken from [5], corresponding 
to a SU(3) quark-gluon plasma with 2+1 light quark flavors, 
have been added for comparison (circles). 



5.3 Adjoint quarks and SUSY Yang-Mills 

Among the various gauge theories built on the SU(A^c) 
gauge group, a particularly interesting case is the one with 
one flavor of massless fermions in the adjoint representa- 
tion. Indeed, such a theory is supersymmetric and it is 
nothing else than a. N = I SUSY Yang- Mills theory [50] . 
the adjoint quarks being also called the gluinos. Within 
our framework, the EoS of that supersymmetric theory 
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can be computed by taking the following pressure as an 
input: 

PSUSY ^Pg+Pq{0), (15) 

with the quarks in the adjoint representation and /i ~ 
0. The observation of the equations presented in Sec. [5] 
reveals that the fermionic contributions to the thermal 
masses only depend on the factors n f and n f dimp re- 
spectively. They are equal to Nc and {N^ — 1) respectively 
in the present case and lead, as in the pure glue case, to 
an EoS of the form {N^ - I) f{T). The ratio A/psB is 
thus predicted to be independent of Nc for A/" = 1 SUSY 
Yang-Mills. 

Figure O clearly shows that the SUSY Yang-Mills EoS 
is identical to the one of the large- iVc QCDas hmit with 
Nf — Uf = 2. Again, the observation of our formalism 
helps to understand that feature, since 



rp I 7V„ — 2 



diiJinnflg^j^Y = - 1, dimRn/l^^^^^ = ^"^^^ ^^ rif. 

(16) 

All those factors are equal at large Nc if Uf = 2, leading to 
identical thermodynamical properties. This provides an- 
other explicit check of the equivalence shown in [4 7) . Two 
remarks remain to be done. First, at large N^ a theory 
with 2 light quarks in the symmetric representation | | | 
(such as used in technicolor models [49|) would have the 
same EoS as SUSY Yang-Mills. Second, to our knowledge, 
the lattice EoS of the TV = 1 SUSY Yang-MiUs theory re- 
mains to be computed. 



5.4 Comments on heavy mesons 

Not only free heavy quarks are of physical relevance to de- 
scribe the QGP, but also heavy mesons. It is indeed widely 
believed that the color interactions are strong enough at 
T > Tc to bind the heavy quarks into mesons. More in- 
formation can be found for example in the review |51|. 
For completeness, let us make some comments about the 
possible large- A^c generalisations of heavy meson models. 

The static qq potential is the basic ingredient of a stan- 
dard heavy meson model. It can be either the qq free en- 
ergy, Fi, or internal energy, Ui = Fi — TOtFi (see e.g. [52] 
for a discussion concerning that choice). Assuming as in 
Sec.[3]that the free energy above Tc is mainly of one-gluon- 
exchange form, as suggested by recent analytical pertur- 
bative calculations |53) and by lattice computations [M,, 
one has -Fi cx —CnX/Nc for a qq pair in a color singlet. 
Comparing the Ci^-factors with the fundamental and A2 
representations, one finds 



a 



N^ 



2N 



1 {Nc + l){Nc~2) 



Nc 



(17) 



Those factors are equal if A^c = 3, but at large- A'c they 
differ by a factor 2, leading to a much more attractive 
potential in the QCDas hmit. In the case of heavy mesons 
also, this last limit is thus expected to exhibit a behavior 
quite far from QCD. 



6 Conclusions 

We have proposed a minimal quasiparticle approach to 
describe the deconfined regime of a generic Q CD-like the- 
ory, with SU(A''c) gauge group and quarks in an arbi- 
trary representation. The model is minimal in the sense 
that we have chosen to use the thermal quark and gluon 
masses computed by using perturbative methods in finite- 
temperature QCD, without extra modifications. The pa- 
rameters are the quark masses, taken from the PDG, the 
critical temperature, for which a standard value is used, 
and the yl-parameter of the running coupling constant. It 
can be fitted on the pure glue case and then taken as an 
input in all the other cases. Its optimal value is found to 
agree with the naive expectation A/Tc = 0{1). 

Despite our model's simplicity, we have first shown 
that it is able to satisfactorily reproduce the equations 
of state computed in lattice QCD well above Tc (typi- 
cally above 2 Tc) either in the pure glue case, or with 
light quark flavors at zero chemical potential. When heavy 
quark flavors or a nonzero chemical potential are taken 
into account, no lattice data are known at high tempera- 
tures. However, since our model may be working well in 
those cases also, it is worth mentioning that it seems to 
extrapolate coherently the existing lattice data. We hope 
that lattice data at nonzero chemical potential and higher 
temperatures will be soon available in order to check our 
model's predictions. 

Then, we have particularized our framework to the 
case of three inequivalent large- A^c limits of QCD: 't Hooft, 
Veneziano, and QCDas- It can be observed that, in the 
pure glue sector, all those limits acts the same since in that 
case no quark flavors are present and the thermodynami- 
cal quantities normalized to the Stefann-Botlzmann pres- 
sure are explicitly independent of Nc at the dominant or- 
der. The universality of the results with respect to Nc has 
been observed on the lattice i4j for A'c > 3 up to the error 
bars. That universality seems to hold even for A'c = 2 at 
T >2 Tc, but not at lower temperatures. Since the SU(2) 
EoS is more than twenty years old [18^, it might be of in- 
terest to make new computations at Nc = 2 with larger 
lattices to see whether the observed low-temperature dif- 
ference is physical (and thus a consequence of the correc- 
tions due to gluon-gluon interactions) or a lattice artefact. 

When light quark flavors are added, the three consid- 
ered large- A^c hmits are inequivalent. It appears that the 
Veneziano limit, followed by the 't Hooft limit, reproduces 
at best the QCD equation of state with light quarks and 
zero chemical potential. It also appears that the change 
in thermodynamical properties due to a nonzero chemical 
potential is an observable which is particularly sensitive 
to the considered large-A'c limit. The QCDas limit, thus 
when quarks are in the two-index antisymmetric represen- 
tation, has the same qualitative behavior as QCD when 
only two light quark flavors are present, but a poor quan- 
titative agreement. With three light quark flavors how- 
ever, even the qualitative agreement is lost. We think that 
this change of qualitative behavior occurs because the cor- 
responding gauge theory develops a Banks-Zaks infrared 
flxed point, absent in QCD |49) . 
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Finally, we have explicitly checked that the equation of 
state in the QCDas hmit with two massless quark flavors 
is equal to the one of the N —1 SUSY Yang-Mills theory, 
i. e. one massless quark flavor in the adjoint representation. 
That last supersymmetric theory begins to be studied on 
the lattice at T = [5S] ; we hope that calculations at T > 
will soon be available to confirm our model's predictions. 

It is worth mentioning that such a quasiparticle pic- 
ture should not, in our opinion, be applied at lower tem- 
peratures unless color interactions are explicitly added to 
the formalism. The maximal value of Z\/T*, systemati- 
cally occurring at T « 1.2 Tc, is a clear indication that 
an ideal-gas approach is questionable near Tc and that an 
interacting quantum gas formalism should be used. We 
leave that issue for future works. 
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